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Abstract 

This paper considers the allocation of time slots in a frame, as well as power and rate to multiple 
receivers on an energy harvesting downlink. Energy arrival times that will occur within the frame are 
known at the beginning of the frame. The goal is to optimize throughput in a proportionally fair way, 
taking into account the inherent differences of channel quality among users. Analysis of structural 
characteristics of the problem reveals that it can be formulated as a biconvex optimization problem, and 
that it has multiple optima. Due to the biconvex nature of the problem, a Block Coordinate Descent 
(BCD) based optimization algorithm that converges to an optimal solution is presented. Numerical and 
simulation results show that the power-time allocations found by BCD achieve a good balance between 
total throughput and fairness. 
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I. Introduction 

Management of energy consumption is vital for the sustainability of many wireless commu- 
nication systems. Especially in the past decade, with the advent of ad-hoc wireless networks, 
energy efficient scheduling policies have been investigated 0], GD> 0. Due to recent advances in 
the areas of solar, piezoelectric and thermal energy harvesting, emerging communication devices 
have been powered by rechargeable batteries which are capable of harvesting energy through solar 
cells, vibration absorption devices, thermoelectric generators, wind power, etc. Although energy 
harvesting allows sustainable and environmentally friendly deployment of wireless networks, it 
requires efficient utilization of time-varying energy. Hence, it shifts the nature of energy-aware 
solutions from minimizing energy expenditure to optimizing it over time. 

It is well known (e.g., flU, (5), (61) that optimization of a broadcast channel (e.g., the downlink) 
shared by many users calls for different choices of rate and power allocation to different users 
depending on the gains, channel conditions, demands of these users, and most importantly, the 
objective of the optimization. We pose the following problem whose objective is proportional 
fairness among users: How to allocate among users the transmission power and the proportion 
of the time between energy harvests, to achieve a good balance between throughput and fairness 
in an energy harvesting broadcast system. Specifically, we investigate the proportional fairness 
based utility maximization problem in a time-sharing multi-user additive white Gaussian noise 
(AWGN) broadcast channel, where the transmitter's battery gets recharged periodically (at known 
intervals). The transmitter wishes to be fair in delivering data to each receiver during the 
scheduling period. Energy is assumed to be harvested at the transmitter during the course of 
transmission. The data, on the other hand, is assumed to be ready at the transmitter before 
the transmission starts. In general, the arrival processes for the data and the harvested energy 
can be formulated as stochastic processes [|7], fl8], requiring an on-line solution that adapts 
transmission power and rate in real-time. However, as discussed by Yang and Ulukus in dU, 
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this is less tractable. We focus on finding the optimum off-line schedule, by assuming that the 
energy arrival profile at the transmitter is deterministic and known ahead of time in an off-line 
manner for a time window, called frame, i.e., the energy harvesting times and the corresponding 
harvested energy amounts are known at the beginning of each frame. The most challenging 
aspect of this optimization problem is the set of causality constraints introduced by the energy 
arrival times, i.e., energy may not be used before it is harvested. 

There has been considerable recent research effort on optimizing data transmission with an 
energy harvesting transmitter. In [TTOl , the authors develop a packet scheduling scheme that 
minimizes the time by which the energy harvesting transmitter delivers all packets to the receiver 
of a single-user communication system. In [11J, the authors extend this work to the multi-user 
case and, propose an iterative approach that reduces the two-user broadcast problem into a single- 
user problem as much as possible, and then, utilizes the single-user solution in IfTOl . f[T2l treats 
the time minimization problem for the two-user broadcast channel differently, as it proposes an 
iterative solution technique by considering two energy arrival slots at a time. These approaches 
are extended by ITT3TI and lfT4l to the case of a transmitter with a finite capacity battery. [|~T5l 
extends ifTOll one step further to propose the directional water-filling algorithm, which is able to 
find the optimal energy management schemes for energy harvesting systems operating in fading 
channels, with finite capacity rechargeable batteries. Both lfl4l and [fT5l investigate the following 
dual offline problems; maximizing the number of bits transmitted with a given deadline constraint, 
and minimizing the transmission completion time with a given number of bits to transmit. 

Unlike the broadcast related studies mentioned above, fl9) investigates the dual problems in a 
multiple access communication system. By using the generalized iterative backward waterfilling 
algorithm proposed in [|9), the authors can simplify the transmission completion time minimiza- 
tion problem into convex optimization problems, and solve the overall problem efficiently. lfT6l 
solves the short-term throughput maximization problem for a battery-limited energy harvesting 
transmitter in a single link topology. 
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In 0/71, the authors consider the problem of energy allocation over a finite horizon for point to- 
point wireless communications, taking into account a time varying channel and energy source, 
so as to maximize the throughput. They use dynamic programming and convex optimization 
techniques to obtain the throughput-optimal energy allocation. In 0, Gatzianas et. al. consider 
an infinite-horizon online throughput maximization problem for a rechargeable sensor network. 
The authors propose a queue stabilizing transmission policy with decoupled admission control 
and energy allocation to maximize a function of the long term rate achieved per link. Chen 
et. al. |fT8l claim that infinite-horizon based solutions can be highly inefficient, especially in 
the context of networks with energy replenishment. Hence, unlike 0, [fT8l investigates the 
finite-horizon throughput maximization problem for a rechargeable sensor network. The main 
difference between two work is that, in [18] inefficiencies cannot be made to vanish, and thus, 
new control techniques are needed to mitigate these inefficiencies. 

This work differs from the previously mentioned studies particularly in its aim to maximize the 
throughput in a proportionally fair way, taking into account the inherent differences of channel 
quality among users. The optimal solution of the problem presented in this paper entails not 
only the optimal power/rate allocation to users between energy arrivals, but also the optimal 
time allocation that will maximize the proportionally fair utility function. Due to characteristics 
of the utility function, the problem presented in this paper can be formulated as a biconvex 
problem and thus it has multiple optima. However, this does not prevent the proposed BCD 
algorithm to converge to one of the many optima, and thus, achieve optimal utility. 

We start by describing the system model in the next section. In Section [nil we make the 
problem statement precise, study the mathematical structure of this problem. The BCD algorithm 
is described in Section [IV] followed by a detailed analysis and discussion of the nature of the 
solution found by BCD. In Section |Vl we test the insight gained from analysis about convergence 
and the nature of the solution, by running the algorithm on numerical examples. We conclude 
in Section [VI] with an outline of further directions. 
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II. System Model 

Consider a time-slotted system where each frame, of length Fi, is divided into K slots. There 
is a single transmitter that transmits to N users by time sharing. The bandwidth for a single link 
channel is W whereas power spectral density of the background noise is iV . Channel conditions 
remain constant during the frame (i.e., g n , the gain of user n, is chosen to be constant throughout 
the frame). The transmitter is equipped with a rechargeable battery such that some energy, E u , 
is harvested from the environment at the beginning of each time slot t of frame i. The length 
of the t th slot of frame i is T ti . In the rest, as we are interested in a specific frame, we drop the 
frame indicator i and define the harvested energy in slot t as E t , and, the length of slot as T t , 
as shown in Fig. [TJ Note that, the slot lengths do not necessarily need to be equal as the energy 
arrivals may occur in different moments in time. We do not restrict our problem formulation to 
the case of periodic energy arrivals ( T t = T for all t G {1, . . . , K} ). In Ifl9ll however, we use 
periodic energy arrivals assumption to derive the characteristics of the optimal solution of the 
problem described in this paper, Problem [1] In this offline problem, energy arrival times and 
amounts that will occur within the frame are known at the beginning of the frame. For a given 
frame, the transmitter chooses a power level p t and a time allocation vector r t = (r u , ...,Tm), 
for each time slot t of the frame, where p nt = p t is the selected transmission power for user n 
during slot t and, r nt is the time allocated for transmission to user n during slot t. 

III. Problem Statement and Structure 

We define the total achievable rate for user n (the total number of bits transmitted to user n 
within the i t h frame), R n = J2?=i T ntW log 2 ^1 + j . Our goal is to maximize a total utility, 
selected as the log-sum of the user rates J2n=i l°92{Rn), which is known to result in proportional 
fairness [20J. However, due to the nature of the time and power allocation problem, and, energy 
harvesting procedure, some constraints need to be satisfied when maximizing the utility function. 
Accordingly, we define the following constrained optimization problem, Problem \T\ where (Q3 
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represent the nonnegativity constraints. The set of equations in ©, called time constraints, ensure 
that the total time allocated to users does not exceed the slot length. The set of equations in ©, 
on the other hand, are technical constraints included to ensure that every user gets a non-zero 
time allocation during the frame. Finally, the set of equations in ©, called energy causality 
constraints, ensure no energy is transmitted before becoming available. 



Problem 1: 

Maximize: U(r,p) = log 2 V] T„tW \og 2 I 1 



N / K 

9nPt 



n=l \t=l 



N W 



subject to: r nt > , p t > (1) 

N 



r «* = T * ( 2 ) 

71=1 

K 

5> nt >e (3) 
t=i 

t t 

Y,PiTi<J2Ei (4) 



where t = 1,...,K and n = 1,...,N. Please note that, Problem [TJ can be written as a 
minimization problem in which the function to be minimized is —U(r,p). One might hope 
that this problem has a unique minimum. Unfortunately, © is a nonlinear non-convex problem 
with potentially multiple local minima, some of which are also globally optimum. Thus, we can 
only expect that by proper choice of the initial value, our algorithm converges to a stationary 
point that is nearby the true optimum. In order to develop such an algorithm, we first decompose 
the problem into two parts (power allocation, time allocation) and determine some characteristics 
that will be useful in understanding the problem structure better. Fortunately, these characteristics 
lead us to Corollary [TJ which establishes that this is a biconvex problem, which we then exploit 
to determine the most appropriate algorithm for solving it. 
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A. Structure of the Optimal Power Allocation Problem 

In this section, we assume that the time allocation, r, is determined, and try to characterize 
the structure of the optimal power allocation problem for this r. When the only variables are 
power variables, Problem Q] reduces to the following constrained optimization problem: 



Problem 2: 


N 

Maximize: U(p) = ^ f n (p) 

n=l 






t t 

subject to: p t > , '^p l T l < ^ Ej 

i=l i=l 


(5) 



where t = 1, K and, f n is a function of the total number of bits sent to user n: 

f n {p) = log 2 ^2 T ntRnt (6) 



,i=l 



and R nt represents the rate of link n in the t slot: 

R nt = Wlog 2 (1 + L n p t ) where L n = (7) 

Lemma \T\ will be useful to get a handle on the characteristics of the problem. Although we 
claim no originality for the results of the lemma, we provide a proof for completeness. 

Lemma 1: i) Let h 1 , . . . , h K be strictly concave functions of pi, . . . ,px respectively, and, 

Ci, ■ ■ ■ ,c K > 0. Then, / = J^f=i c i^i ^ s concave. If one of the q's is positive (> 0), then / 

is strictly concave. 

ii) Increasing concave functions of strictly concave functions are strictly concave. 

Proof: The proof is provided in Appendix |A] ■ 
Theorem 1: Problem [2] can be formulated as a strictly convex optimization problem. Thus, 
there exists only one global optimum for a given time allocation. 

Proof: The proof is provided in Appendix |Bj ■ 
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B. Structure of the Optimal Time Allocation Problem 

In this section, we assume that the power allocation across all slots has been determined. 
Then, given that the power variables are known constants, we determine the characteristics of 
the time allocation. So Problem \T\ reduces to Problem [3] where the only variables are the time 
variables: 



Problem 3: 










N 








Maximize: U (r) = s n (r) 

n=l 








N 


K 






subject to: r nt > , ^ T nt = T t , 


^ T nt > £ 


(8) 




n=l 


t=l 





where t — 1, K, n — 1, N and, s n is a function of the time variables: 



s n (r) = log 2 ^^r nt i? nt J (9) 

and Rnt's (defined in Eq. |7]) are known constants that represent the rate of link n in the t th slot. 

Theorem |2] below records the convexity of Problem 3, and leads to one of the main results of 
this paper, Corollary [TJ The proof of Theorem |2]rests on the observation in Lemma |2l Although 
we claim no originality for the results of the lemma, we provide the details for completeness. 
Lemma 2: i) Let qi, . . . , qx be affine functions of r nl , . . . , r nK respectively, and d\, . . . , &k > 

0. Then, m = J2f =1 is affine. 
ii) Increasing concave functions of affine functions are concave. 

Proof: The proof is provided in Appendix ■ 
Theorem 2: Problem |3] can be formulated as a convex optimization problem. Thus, all local 
optima are global optima. 

Proof: The proof is provided in Appendix |D] ■ 
Note that, Problem |3] is convex, but not necessarily strictly convex. Therefore, in general, rather 
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than a unique global optimum, there may be multiple local optima which are all also globally 
optimum. 

Corollary 1: Problem \T\ can be formulated as a biconvex optimization problem. 

Proof: The proof is provided in Appendix |E] ■ 

IV. Solution Method 

In the previous section, we have shown that Problem Q] can be formulated as a biconvex 
optimization problem since —U(r,p) is a biconvex function. Such functions are well-studied in 
the optimization literature ED . [1221 . While not convex, they admit efficient coordinate descent 
algorithms that solve a convex program at each step. In this section, we present a block coordinate 
descent based algorithm, shortly BCD, for solving Problem [TJ In the BCD solution method, 
sequentially one block of variables is minimized under corresponding constraints while the 
remaining blocks are fixed. We have the simplest case of only two block variables r and p. 
Hence, the algorithm alternates between minimization with respect to r and minimization with 
respect to p. Our BCD algorithm operates explicitly as follows: 

1) Start from any valid time allocation, for example assign each time slot to different user in 
the form of TDMA. Assuming that all of the energy E t is used up until the end of period 
t, the power is determined. This power setting satisfies Eq. ©. 

2) Keep r nt fixed for all n and t. Optimize U(r,p) with respect to p t , t — 1, . . . , K and the 
constraints given by ©. 

3) Repeat the following for all t — 1, . . . , K.; Keep r ni fixed for all n = 1, . . . , N and i ^ t. 
Also keep p t fixed for all t. Maximize U(j,p) with respect to T nt ,n = 1,...,N and 
constraint in Eq. ©. 

4) If the variables have converged, stop. Otherwise, go to Step 2. 

For optimization of the time variables, the Lagrange multiplier method is used. The opti- 
mization of the power variables, however, is accomplished by using the Sequential Uncon- 
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strained Minimization Technique (SUMT) 11231 . SUMT is an optimization method that converts 
a constrained optimization problem into an unconstrained one by adding the constraints to the 
objective function as a "penalty". It then uses a standard unconstrained optimization algorithm 
(e.g., Newton, Steepest Descent, etc.) Il24l . Il25l to solve the problem with the new objective 
function. 

Regarding the issue of convergence, Problem [1] is a biconvex optimization problem and as 
such potentially, there exist many local optima. Therefore, convergence to the global optimum 
is not guaranteed. However, provided that some conditions are satisfied, convergence to a partial 
optimum (see Definition Q]) is guaranteed. As also discussed by Lin in G6l . convergence to a 
stationary (or critical) point, for block coordinate descent methods requires sub-problems to have 
unique solutions ( 11271 , [|28l ), but this property does not hold here: Although sub-problem |2] is 
strictly convex, |3] is not strictly convex (only convex). Fortunately, for the case of two blocks, 
Grippo and Sciandrone [|29l have shown that this uniqueness condition is not needed. Hence, 
BCD converges to a stationary point of Problem [TJ As a stationary point can be; minimum, 
maximum, or a saddle point, this convergence result may not be sufficient. However, we can still 
use the following definition and theorem (Definition 4.1 and Theorem 4.2 of [1271 . respectively) to 
build a stronger result. For this, let X C 9ft n and Y C !R m be two nonempty sets, let B C X x Y, 
and, let B x and B y denote the x-sections and y- sections of B, respectively. 

Definition 1: Let / : B — > 3? be a given function and let (x*, y*) G B. Then, (x*, y*) is called 
a partial optimum of / on B, if 

f(x*,y*)<f(x,y*) \/xeB* y and f(x*,y*)<f(x*,y) Vy G B* x (10) 

Theorem 3: Let B be a biconvex set and let / : B — V 3ft, be a differentiable, biconvex function. 
Then, each stationary point of / is a partial optimum. 
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Furthermore, Theorem 4.9 of [21] shows that, for BCD-like algorithm!] (There are only two 
block of variables, and, sequentially one block of variables is minimized under corresponding 
constraints and the other block is fixed), when one of the subproblems is strictly convex, all 
partial optima have the same function value. Hence, we conclude that the BCD algorithm surely 
converges to a partial optimum of Problem [H and all partial optima of Problem \T\ yield the same 
utility value. Note that although the final allocation, (r*,p*) generated by the BCD algorithm, 
might be a partial optimum, it neither has to be a global nor a local optimum to the given 
biconvex optimization problem. The following theorem by Wendell and Hurter [|30l describes the 
connection between partial and local optima. For the following biconvex minimization problem, 

min{f(x,y) : x G X C W 1 , y G Y C 3? m } (11) 

Theorem 4: Let X C 9ft n and Y C be convex sets and let / : B — > be a biconvex 
function with a partial optimum in (x*,y*) G X x Y. Also, let U(y*) be the set of all optimal 
solutions to the problem in (fTT|) with y — y* and let U(x*) denote the set of all optimal solutions 
to (PTTI) with x = x*. If (x*, y*) is a local optimum of (fTTT) . then it necessarily holds that 

f{x*,y*)<f{x,y) VxeU(x*) VyeU(y*) (12) 

We have observed through extensive numerical analysis that the above condition is invariably 
satisfied for the partial optimum found by BCD, which strengthens our belief that BCD achieves 
optimal or close to optimal utility. However, as also noted in fl2TJ, the given local optimality 
condition is in general not sufficient. Indeed, Wiesemann claims in ||3TI (p. 92) that, even the 
verification whether a particular solution to a biconvex problem is locally optimal is NP- 
complete. Gorski et.al. [21], on the other hand, claims that to find the global optimum of a 
biconvex minimization problem by a BCD-like algorithm (ACS ETTl ). a multistart version of 

'ACS (Alternate Convex Search) algorithm proposed in 11211 . 
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BCD can be used. But, still, there is no guarantee to find the global optimum within a reasonable 
amount of time or to be sure that the actual best minimum is the global one. Hence, it seems 
justified to settle for the modest goal to find a partial optimum in our case. 

V. Numerical and Simulation Results 

In this section, we present the numerical and simulation results related to BCD algorithm. 
Throughout our simulations we use the folowing setup: W = 1kHz, N Q = 10~ 6 W/Hz. 
Unless otherwise stated, all powers are in Watts and all energies are in Joules. For the sake 
of an example, we suppose that there are five users in the system and 10 energy arrivals in 
100 sees (frame length). The arrivals are E = [20,100,1,1,1,70,100,1,10,40] joules in the 
[I s *, 2 nd , . . . , 10 th ] slots respectively. The first user is the strongest one, and, other users are 
ordered in a such way that the preceding user is twice as strong as the previous one, i.e., path 
losses of the users are; 25, 28, 31, 34, 37 dB respectively. The starting point of the algorithm 
is the "Spend What You Get" policy (proposed by Gorlatova et. al. (H) combined with TDMA 
time allocation. This policy corresponds to using all energy in the epoch it was harvested in, and 
will be referred to in the rest as SG+TDMA. We performed simulations both for unequal and 
equal slot lengths. In our simulations, we use the following sequence of slot sizes for 10 slots; 
Si = [10, 12, 5, 7, 4, 15, 20, 2, 10, 15] and S 2 = [25, 44, 14, 7, 3, 32, 47, 19, 26, 38], for the case of 
unequal slot lengths, and, Si = [10, 10, ... , 10] and S 2 = [25.5, 25.5, . . . , 25.5] for the case of 
equal slot lengths. Note that, Si and Si have the frame length of 100 sees, whereas S 2 and S 2 
have the frame length of 255 sees. 

First, we assume that the frame length is 100 sees, and, we illustrate the power iterations of 
the BCD algorithm, for Si, in Figure |2] The power convergence of the algorithm for periodic 
energy arrivals (Si), however, is illustrated in Figure |2] As observed from the figures, rather 
than transmitting with full power, saving some energy for the future use is preferred. Another 
observation is about the fast convergence of the algorithm, i.e., the powers seem to rarely change 
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after just a few iterations. In order to observe the effect of frame length and different slot lengths, 
in Figure @] we show how utility improves through the iterations for all of the aforementioned 
slot length sequences, i.e., Si, Si, S2, S2. The fast convergence of the BCD algorithm is more 
evident in this figure. The optimal schedules (power and time), optimal utility and thus, the 
utility improvement (when compared to SG+TDMA) obtained by BCD, for all four sequences 
are presented in Table IB 

In some energy harvesting systems, transmitters have supercapacitors that can store the har- 
vested energy and supply in every predetermined time window, allowing the case of periodic 
energy arrivals. In such a case, if no energy is harvested within a slot, we set the amount of 
harvested energy to for that slot. As observed from Figure Hand Table U periodic energy arrivals 
assumption does not degrade the system performance, yet increases the utility improvement. 
Moreover, as we have shown in [fT9l , by using the periodic energy arrivals assumption we 
can analytically derive the characteristics of the optimal solution of Problem [TJ and, develop two 
heuristics that closely track the performance of BCD algorithm. Hence, from now on, we present 
results only for the case of periodic energy arrivals in 100 sees, Si. 

Throughput improvement is another important criteria in our problem setup. Hence, we next 
investigate the throughput improvement of the users for increasing path losses. The results are 
illustrated in Figure |5] In the figure, the Mean Path Loss, is computed as L = i J2iLi Li 
where Li represents the path loss of user i. As seen from the figure, with minor decrease in 
the throughput improvement of the stronger users, the weak users receive much more bits than 
that they used to receive with SG policy and TDMA. As L increases, the overall throughput 
improvement also increases. For instance, when L = 31, User 1, User 3, User 4, and, User 5 
enjoy approximately 3 %, 1621 %, 361 %, 80 % throughput improvement respectively, while 
User 1 suffers only 32 % of loss. Clearly, BCD is a proportionally fair algorithm which tries to 
maximize the utility by meeting certain demands of every user. 

We next analyze the effect of number and amount of energy arrivals. Assume that there are six 
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users in the system with the following path losses: 19, 22, 25, 38, 31 dB. The results obtained 
for six different energy arrival sequences (two for K=10 and four for K=12) are shown in Table 
HO The arrival sequences are intentionally chosen similar to each other, so that it would be 
easier to determine the effect of small changes on the utility improvement. As observed, it is 
not the number of slots (number of energy arrivals) but the amount of every individual energy 
harvest that determines the utility improvement. The events like; sudden decrease in energy level 
or harvesting very small amount of energy for a long time increases the utility improvement 
obtained by BCD, as SG+TDMA policy may cause the base station to stay idle for a long time 
because it does not save energy for future use. 

Knowing that it is not the number but the nature of harvests that affect the utility improvement, 
from now on, we set the number of energy arrivals to be 10 (K=10) and the harvests to arrive 
as in the 2 nd arrival sequence, i.e. [20, 100, 1, 1, 1, 70, 100, 1, 10, 40] so that we can analyze 
the effect of number of users to the performance of the BCD algorithm. Keeping the number of 
harvests and harvest values the same, we perform a series of simulations with different number of 
users. First, the effects of the optimal power-time allocation pairs on utility, utility improvement, 
and fairness are investigated. In order to be able to analyze all scenarios, in the next three figures, 
we use the following setups: a) The strongest user in the system has 13 dB path loss, and, every 
new user that joins the system has 3 dB more path loss than the previous one. b) The strongest 
user has 19 dB path loss, and, every new user deviates 3 dB. c) The strongest user has 25 dB 
path loss, and, every new user deviates 3 dB. Hence, 13 dB, 19 dB, and, 25 dB seen in the 
figures represent the path loss of the strongest user in the system. Figures [6] and [7] show how 
utility and percentage improvement in utility, respectively, change with the increasing number 
of users. As seen from the figures, the solution found by BCD exhibits significant improvement 
over a SG+TDMA schedule. Between two methods, SG+TDMA is the worst one since even 
with a few users, utility can be improved. The results show that when case c) is valid, a utility 
improvement of approximately 20% is possible with BCD. 
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Although we aim at proportional fairness in this work, it may be interesting to analyse max- 
min fairness of the BCD algorithm. Jain's index [1321 . Il33l is a well-known measure of fairness. 
The index FI takes the value of 1 when there is a complete fair allocation. 

FI = (Eti^) 2 (13) 

For computing FI, we use the no. of bits transmitted to the users, Xi = 2 Ui for i = 1, . . . , N, 
where U t = log 2 (j2%=i T it R it) ■ From Figure E it is clear that SG+TDMA is worse than BCD 
in terms of fairness. Especially for eight users, FI S g+tdma = 0.41 whereas FI BC d = 0.80. 
Although low path losses embrace lower utility improvement, they mainly allow BCD algorithm 
to be very efficient in terms of fairness, e.g., above 0.8. However as the path loss difference 
between users increase, completely fair allocations may not be the optimal ones. For instance, 
when there are eight users in the system the path losses of the users in case c) are; 25, 28, 31, 
34, 37, 40, 43, 46 dB, yielding an excessive difference of 23 dB between the weakest and the 
strongest user. In this case, the algorithm should favor user 1 more than it favors user 8, in order 
to maximize the utility function, causing a proportionally (instead of purely) fair allocation. 



VI. Conclusion 

This paper investigated the proportional fair power and time allocation problem in an energy 
harvesting broadcast system. The paper focuses on finding the optimum off-line schedule, by 
assuming that the energy harvesting times and the corresponding harvested energy amounts are 
known at the beginning of each frame. In order to determine the optimal schedule, the problem 
is designed as a proportional fairness based utility maximization problem. Detailed analysis 
of structural characteristics of the problem has been performed, which revealed that it can be 
formulated as a biconvex optimization problem, and that it has multiple optima. Furthermore, an 
algorithm based on block coordinate descent (BCD), that surely converges to a partial optima 
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of the problem, has been showed. Simulation results show that by allocating among users the 
transmission power and the proportion of the time between energy harvests, BCD achieves a 
good balance between throughput and fairness on an energy harvesting downlink. BCD algorithm 
converges to the partial optimal utility and reaches up to 20% utility improvement compared to 
SG+TDMA schedule. 

Ifl9l will investigate the power and time related characteristics of an optimal solution of the 
proposed problem, and, develop low-complexity heuristic algorithms that will closely track the 
performance of the BCD algorithm. The analysis, also, will be validated by the consistency of 
the numerical results and derivations. 

Appendix A 
Proof of Lemma [j] 

i) Let us define I = J2f=i c i^i where hi is a strictly concave function of Pi, and, > 0. 
For I to be concave, it needs to satisy the concavity condition, i.e., l(Xpi + (1 — X)p^) > 
Xl(pi) + (1 — A)i(pi) where < A < 1, for any point pi, p2 in the domain of I. Thus, 
proving that / satisfies the concavity condition completes the proof of part (i). We start by 

K K 

l{Xpi + (1 - A)pi) = c A(Api* + (1 - X)p 2i ) > J2 c i( Xh i(Pu) + (! - A)/k(p 2 i)) (14) 

i=i i=i 

K K 

= XJ2 Cihiipu) + (1-X)J2 cMP2i) = AZ(pT) + (1 - (15) 
i=i t=i 

where Eq. (fl4l) follows from the strict concavity of hi function. From the set of equations 

described above, one can observe that I is a concave function of p. Note that in Eq. (fl4l) . 

equality may happen only if q = for all i. Hence, if there exists an index j such that 

Q = for i — 1, . . . , j — 1, j + 1, . . . , K and, q > for i — j, Eq. (fT4l) is satisfied with 

">". Then, I is strictly concave in pj:. This completes the proof of part (i). 
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ii) This is rather a general result, and hence, due to space constraints, we refer the interested 
reader to our technical report [34] for the details of our version of the proof. 



Appendix B 



Proof of Theorem [j] 



Since the constraints of the problem are linear, and, maximizing U(p) is equivalent to min- 
imizing —U(p), showing that the utility function, U(p), is strictly concave will be enough to 
show that Problem [2] can be formulated as a strictly convex optimization problem. We start by 
checking the concavity of R nt . As R nt is a function of p t , let h t = R n t (as in the proof of 



negative for all t = 1, . . . , K. Hence, from the second derivative test ||23~1 , h t , thus R nt , is 
strictly concave in p t . Furthermore, let l n = J2tLi T ntRnt for an arbitrary user n. As all r ni 's 
(for t = 1, . . . , K) are nonnegative, and, at least one r nt is positive for user n, from part (i) of 
Lemma [H /„ is a strictly concave function of p. Note that this is true for all users, i.e., l n is 
strictly concave for all n = 1, . . . , N. Thus, from part (ii) of Lemma [H /„ is strictly concave 
in p for all n = 1, . . . , N. The rest of the proof is straight-forward, since the utility function, 
U(p), is a positive linear combination of / n 's and thus (from part (i) of Lemma [D is strictly 
concave in p. Hence, the proof is complete. 



i) Let us define m = J2i=i d-iQi where q; t is an affine function of r ni , and, rfj > 0. For 
m to be affine, it needs to satisfy the affinity condition, i.e., m(\T[ + (1 — A)tJ) = 
Am(rT) + (1 — \)m(T^), where < A < 1 and Tji: = [r fc nl . . . T kynK ] T is a point on m. 
Thus, proving that m satisfies the condition completes the proof of part (i). We start by 





^ is definitely 



Appendix C 



Proof of Lemma [2] 
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K K 

m(Xn + (1 - A)t5) = diqi(\T U + (1 - A)r 2l ) = ^^(Ag^n*) + (1 - A)%(r 2 i)) (16) 

i=l i=l 
K K 

= \J2 d Mi(ni) + (1-X)J2 diUfoi) = Am(7T) + (1 - A)m(^) (17) 

i=l i=l 

where ru, r 2i are the i th entries of the rT vector and vector respectively, and, Eq. (fT6l 
follows from the "affine" property of the rrii function. From Eq. ([171 . one can clearly 
observe that m function is an affine function of time variables. Hence, any nonnegative 
linear combination of affine functions is affine, and, the proof of part (i) is complete, 
ii) Similar to part (ii) of Lemma [T] this is a general result that can be found in the literature. 
Hence, due to space constraints, we refer the interested reader to ll34l for the details. 



Appendix D 
Proof of Theorem [2] 



Similar to the proof of Theorem [T] showing that the utility function, U(r), is concave will 
be enough to show that Problem [3] can be formulated as a convex optimization problem. We 
start by checking the concavity of s n for n = 1, . . . , N. As R nt is a known constant for every 
n E {1, . . . , iV} and t G {1, . . . , K}, let d t = R n t, and, q t = r nt for an arbitrary user n (as in 
the proof of Lemma |2fs part {%)). It is well-known that a linear function is an affine function. 
As q t is a linear function of r nt for any user n, q t is affine. Now, let m n = J2?=i T ntR n t for an 
arbitrary user n. Note that all i? nt 's are nonnegative constants known apriori. Thus, from part (ii) 
of Lemma [2l s n is concave in r for all n = 1, . . . , N. The rest of the proof is straight-forward, 
since the utility function, U(t), is a nonnegative linear combination of s n 's and thus (from part 
(i) of Lemma CO) is concave in r. Hence, the proof is complete. 
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Appendix E 
Proof of Corollary Q] 

A function / : X x Y — > 3ft is called biconvex if f(x, y) is convex in y for fixed x E X and is 
convex in x for fixed y E Y Il22l . Since the constraints of the problem are linear, showing that 
—U (r, p) is biconvex will be enough to show that Problem \T\ can be formulated as a biconvex 
optimization problem. —U(r,p) is a function of two set of variables; r and p. From Theorem 
Q3 given r, —U(p) is convex. Similarly, from Theorem |2] given p, —Uir) is convex. Hence, 
—U(t,p) is biconvex, which completes the proof. 
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Fig. 1. Problem illustration: There are K energy arivals in a frame, and, the time between consecutive arrivals are allocated 
to N users. 
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TABLE I 

The Results of BCD Algorithm for four different slot length sequences. 
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Iterations 

Fig. 2. Powers vs. Iterations (N=5, K=10, Unequal slot lengths): The numbers in the legend represent the corresponding slots. 
Starting from SG policy, BCD converges to the optimal powers in 1 1 iterations. 



TABLE II 

The effect of number and amount of energy harvests 
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Fig. 4. Utility vs. Iterations (N=5, K=10): Starting from SG+TDMA, BCD converges to the optimal utility in 11,8,16,18 
iterations for the following slot length sequences respectively: Si, Si, S2, S2. 
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Fig. 5. Throughput Improvement vs. Mean Path Loss (N=5, K=10): Mean path loss is computed as the mean of the path losses 
of all users in the system. Results represent the throughput improvement of five users for three different path loss patterns. With 
minor decrease in the throughput of the stronger users, the weak users receive much more bits than that they used to receive 
with SG+TDMA. 
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Fig. 6. Utility (SG+TDMA, BCD) vs. No. of Users: The utilities obtained by SG+TDMA and the proposed algorithm, for 
increasing number of users, are compared. The effect of path loss (the strongest user's path loss is shown between parentheses) 
on utility is shown. As path losses of the users increase the utility decreases. 
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Fig. 7. Utility Improvement vs. No. of Users: The utility improvement of the proposed algorithm over SG+TDMA, for increasing 
number of users, are compared. The effect of path loss (the strongest user's path loss is shown between parentheses) on utility 
improvement is shown. As path losses of the users increase the utility improvement increases, i.e, BCD's performance improves 
as the channel quality becomes degraded. 
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Fig. 8. Fairness Index (SG+TDMA, BCD) vs. No. of Users: The fairness of SG+TDMA and the proposed algorithm, for 
increasing number of users, are compared through FI, which takes the value of 1 when there is a complete fair allocation. The 
effect of path loss (the strongest user's path loss is shown between parentheses) on fairness is shown. As difference among the 
path losses of the users increase the fairness indexes of the schemes decrease (SG+TDMA being the worst), causing comparatively 
unfair allocations among users. 
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